Spectroscopic LSJ notation for atomic levels obtained 

Q ; from relativistic calculations 

o 

(N 

^ ! G. Gaigalas "'''* T. Zalandauskas '^ and S. Fritzsche" 

1—5 

^r^ \ °- Fachbereich Physik, Universitat Kassel, 

Heinrich-Plett-Str. 40, D-34132 Kassel, Germany. 
^', 

^ Vilnius University Research Institute of Theoretical Physics and Astronomy, 
A. Gostauto 12, Vilnius 2600, Lithuania. 

^ : February 2, 2008 

O 

•i-H 

>-. 

(~| ■ Abstract 

Today, relativistic calculations are known to provide a very successful means in the 
study of open-shell atoms and ions. But although accurate atomic data are obtained from 
^ , these computations, they are traditionally carried out in jj— coupling and, hence, do often 

^«0 ' not allow for a simple LSJ classification of the atomic levels as needed by experiment. In 

fact, this lack of providing a proper spectroscopic notation from relativistic structure cal- 



o 

o . . 

\^ , culations has recently hampered not only the spectroscopy of medium and heavy elements, 

^^ ' but also the interpretation and analysis of inner-shell processes, for which the occurrence 

of additional vacancies usually leads to a very detailed fine structure. 

Therefore, in order to facilitate the classification of atomic levels from such computa- 
fj ' tions, here we present a program (within the Ratip environment) which help transform 

'^ , the atomic wave functions from jj— coupled multiconfiguration Dirac-Fock computations 

^->| into a LS"— coupled representation. Beside of a proper LSJ assignment to the atomic lev- 

els, the program also supports the full transformation of the wave functions if required for 
(nonrelativistic) computations. 
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PROGRAM SUMMARY 

Title of program: LSJ. 

Catalogue number: ADTL 

Program summary URL: http://cpc.cs.qub.ac.uk/suinmaries/ADTLl 



Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland. 

Licensing provisions: None. 

Computer for which the new version has been tested: IBM RS 6000, PC Pentium III. 
Installations: University of Kassel (Germany). 
Operating systems: IBM AIX 4.1.2+, Linux 7.1.+ 

Program language used in the new version: ANSI standard Fortran 90/95. 

Memory required to execute with typical data: Memory requirements depend on the shell 
structure and the size of the wave function expansion which is used to represent the atomic 
levels. 

No. of bits in a word: All real variables are parametrized by a selected kind parameter 
and, thus, can easily be adapted to any required precision as supported by the compiler. 
Presently, the kind parameter is set to double precision (two 32-bit words) in the module 
rabs_constant. 

CPU time required to execute test data: 5 seconds on a 1 GHz Pentium III processor. 

Distribution format: compressed tar file. 

Keywords: atomic, LSJ spectroscopic notation, LS — jj transformation, multiconfiguration 
Dirac-Fock, recoupling of angular momenta, relativistic. 

Nature of the physical problem: 

The spectroscopic LSJ notation is determined for atomic levels which were calculated pre- 
viously in the framework of the jj— coupled multiconfiguration Dirac-Fock (MCDF) model. 
This notation is based on a complete jj — LS transformation of the leading jj— coupled con- 
figuration state functions (CSF) in the wave function representation of the selected levels. 

Restrictions onto the complexity of the problem: 

The jj — > LS transformation of the (jj— coupled) CSF is supported for all shell structures 
including open s—, p—, d—, and /—shells. For shells with / > 3 (i.e. beyond the /— subshells), 
however, a proper transformation of the antisymmetrized subshell states can be carried out 
only for the case of one or two equivalent electrons. This restriction also applies for the 
transformation of the 57/2 and gg/2 subshell states which are otherwise supported by the 
Ratip package [S. Fritzsche, J. Elec. Spec. Rel. Phen. 114-116 (2001) 1155.]. The jj ^ LS 
transformation matrices, which are applied internally by the program, are consistent with the 
definition of the (reduced) coefficients of fractional parentage [G. Gaigalas, Z. Rudzikas and 



C. Froese Fischer, At. Data Nucl. Data Tables, 70 (1998) 1; G. Gaigalas, S. Fritzsche and 
Z. Rudzikas, At. Data Nucl. Data Tables, 76 (2000) 235.] as published previously. 

Unusual features of the program: 

The LSJ program is designed as a part of the Ratip package [S. Fritzsche, J. Elec. Spec. 
Rel. Phen. 114—116 (2001) 1155.] for the computation of (relativistic) atomic transition 
and ionization properties. This (new) component therefore supports the transformation of all 
atomic states which are generated either with Ratip or by means of the Grasp92 code [F. A. 
Parpia, C. F. Fischer, and I. P. Grant, Comput. Phys. Commun. 94 (1996) 249.]. Moreover, 
the normalization of the transformed states is tested within the LS'J— coupled basis. 



LONG WRITE-UP 



1 Introduction 

In atomic spectroscopy, the standard LSJ notation of the atomic levels is frequently applied 
for classifying the low-lying level structures of (most) atoms and ions. The knowledge of the 
basic LSJ symmetries of the levels often helps the interpretation of atomic spectra without 
detailed computations. In fact, the LSJ spectroscopic notation is found useful also in order 
to explain the excitation and the decay properties of individual levels, at least qualitatively. 
Since, in addition, most atomic structure calculations were carried out within a LS*— coupled 
basis in the past, the LSJ designation of the levels is taken today as the standard notation in 
the literature as well as in the most data bases. 

In recent years, however, relativistic structure computations also became feasible and were 
frequently performed, when accurate data were required. Based on either the Dirac-Coulomb 
or Dirac-Coulomb-Breit Hamiltonian, for instance, these computations are traditionally car- 
ried out in jj— coupling, from which the LSJ notation cannot be derived so easily, especially 
if atomic configurations with (several) open shells occur. In the last few years, therefore, the 
lack of a proper spectroscopic notation deriving from relativistic calculations has hampered not 
only further analysis of open d— and /—shell elements but also the study and interpretation 
of the inner-shell processes where the occurrence of (additional) inner-shell holes often result 
in a very detailed fine-structure. In order to obtain the standard notation in such a case, an 
additional transformation of the atomic states (from the j'j— coupled into a LS"- coupled basis) 
is needed first before the full benefit is drawn from such a relativistic treatment. But although 
several tools and tabulations are available today for carrying out the transformation between 
different coupling schemes, there is no program available which would directly support the 
well-known relativistic structure codes such as Grasp92 1 or Ratip |51I3|^. The reasons 
for this lack are simple: (i) the rapidly increasing complexity of the transformation matrices, 
when several open shells are involved, and (ii) no simple access to the standard quantities from 
the atomic shell model. 

Apart from the explicit jj — LS transformation of the many-electron basis, there is a second 
method known from the literature 01^1 d which helps to derive the LSJ classification of atomic 
levels from relativistic computations. This method applies the (effective) relativistic operators 
within an L5— coupled basis and, hence, can be directly implemented into the non-relativistic 
theory, i.e. by using the familiar Pauli spinors for the construction of the many-electron basis. 
So far, however, these effective operators have neither been derived in detail for all the physical 
cases of interest jHj nor does a full implementation exist which is competitive with the standard 
relativistic codes. Moreover, the use of effective operators within a Pauli spinor basis cannot 
account for the different radial behaviour of equivalent electrons with j = I ± 1/2 {I > 1) as 



is found to be important, in particular, for the behaviour of inner-shell electrons in medium 
and heavy elements. Therefore, the explicit jj — LS transformation of the many-electron 
basis seems to us a much simpler route to derive the LSJ designation of the atomic states, 
making use of the well-known relativistic structure codes such as Grasp92 pQ or Ratip j^J 
as its input. This jj — LS transformation of the basis states has recently been generalized 
also to allow an arbitrary number of open shells in the construction of the symmetry-adapted 
functions 0. 

To facilitate the classification of the atomic levels and wave functions from the relativistic 
calculations, here we present the Lsj program as a new component of the Ratip package 
[2]. This component supports both the LSJ classification of the level structure of open-shell 
atoms and ions as well as the (complete) transformation of the corresponding wave functions. 
But the program may also be of value in solving the problems of a more general character, 
where the transformation from the jj- into LS'-coupling scheme and vice versa is performed, 
e.g. a transformation of this type is needed when one seeks the optimal coupling scheme after 
calculating the energy spectra in a non-relativistic or relativistic approximation. But the 
problems of this type are quite independent of those investigated here. They are related only 
to a similar transformation being made. Therefore the authors have not aimed at creating a 
general program in the frame of the Ratip package that would be of use while the physical 
problems are formulated in the most general case, especially with certain risks originating 
from such a formulation of a problem when the transformation considered is based on the 
transformation coefficients jJOl- The values of these depend on the phase convention that 
is chosen to obtain the spin-angular coefficients and on the definition of the coefficients of 
fractional parentage (CFP). Of course, the program presented in this work may serve as a 
basis for solving the other problems in atomic physics. But the authors have not aimed at 
the exhaustive description of the program on purpose, because in that case the technical way 
of answering those questions is not as important as is the methodological one, related to the 
phase convention and the definition of CFP. The program is based on the methodology that 
uses the Fano-Racah convention, and the CFPs are defined via the quasispin formalism |llj . 
The CFPs are published in |121 I13j . In all other cases this transformation leads to erroneous 
results. Therefore, before starting to perform any tasks related to the use of the program in 
other packages or trying to solve any problems of other type, we recommend that the authors 
of the paper should be contacted. 

In the following section, we first provide a brief outline on the transformation of the atomic 
states in going from the jj— coupled to a LS"— coupled many-electron basis. This shows how 
the LSJ spectroscopic notation is obtained as a symmetry of the leading term in the wave 
function representation within the LS"— coupled basis; however, we leave out most of the details 
concerning the construction and recoupling of the symmetry functions in the atomic shell model 
since that material has been presented recently ^J|. In section 3, we describe the program 
organization and how the code is distributed. To demonstrate the interactive use of the Lsj 
program, section 4 displays — by two examples — the derivation of the LSJ notation for the 
low-lying excited levels of beryllium as well as for the ground-state of gadolinium. Finally, 
the conclusions and a brief outlook are given in section 5. 



2 Theoretical background 

The task to derive the LSJ designation of atomic levels from relativistic calculations can be 
summarized rather easily without much knowledge about the atomic shell model. Having, 
for example, a representation of an atomic state ^ with respect to a jj— coupled symmetry 
basis, the (first few) leading terms in the wave function expansions need to be transformed 
into a -LS'( J)— coupled basis from which, then, the proper spectroscopic notation is given by 
the symmetry of (again) the leading term. At the same time, the contributions from the 
subsequent terms in the L5— coupled basis indicate the admixture of other LS symmetries 
and, hence, the deviation from a pure LS or Russell-Saunders coupling scheme. Of course, 
the explicit jj — LS transformation of the atomic state functions (ASF) requires quite detailed 
information about the construction and the coupling of the symmetry adapted basis functions, 
both in jj— as well as -L5— coupling. For most practical applications, however, these details 
are irrelevant. 

In this write-up of the Lsj program, we therefore restrict the discussion mainly to the repre- 
sentation of the atomic states in the various coupling schemes, but otherwise refer the reader to 
the literature ^lEl- While this restriction enables us to explain the basis steps of a jj — LS 
transformation of atomic states from standard relativistic computations, and how these steps 
are implemented within the program, it avoids the need of providing too many details about 
the definition of symmetry-adapted functions and about the theory of transformation in the 
atomic shell model. 

2.1 Representation of atomic states 

In principle, of course, any valid coupling scheme can be used in order to represent the (many- 
electron) wave functions in atomic structure computations. In relativistic calculations, how- 
ever, a jj— coupled basis appears particularly appropriate as the symmetry of the Dirac equa- 
tion for a single electron leads naturally to solutions which have a well-defined total angular 
momentum, j = 1 + s, but no defined projections (onto the z— axis) of either the orbital angu- 
lar momentum 1 nor the spin s of the electron. In the multiconfiguration Dirac-Fock (MCDF) 
or configuration interaction (CI) methods, these one-electron solutions of the Dirac equation 
(which are finally determined for some effective potential) are then used also for the set-up of 
the many-electron basis functions, the so-called configuration state functions (CSF). Owing 
to the overall symmetry of atoms with respect to a rotation or inversion of the coordinates, 
these CSF are constructed in order to have a well-defined total angular momentum J and 
total parity P. For a given set of such configuration states, W"frJP) , f = ^■■■nc}, the atomic 
states are obtained simply as a linear combination 

I'^riJ'')) = X^ ato)(r) l7.JP) (1) 

r 

of CSF of the same symmetry, where r = 1, 2, ... enumerates the atomic levels (of the given 
symmetry) and jr denote the set of all further quantum numbers as required for an unique 
classification of the symmetry-adapted basis. Obviously, these quantum numbers and the 
details of their coupling play an important role in the transformation of the basis from one 
coupling scheme into another. For the present discussion, however, it is sufficient to say that 
the CSF basis in (^ spans a (finite) part of the A^— electron Hilbert space and that the Ur {t) 
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represent the (usual) Fourier coefficients for the expansion of an atomic state with respect to 
this basis. As seen from ansatz (pQ), moreover, only J and P are in fact good quantum numbers 
which, apart from the subshell angular momenta (see below) and the total energy of the ASF, 
can be used for the classification of the atomic levels. 

An (almost) equivalent^ but alternative representation of the ASF can be obtained within a 
LS*— coupled basis 

l^-('^'')> = Y.''s'''\^)\lsLSJP), (2) 

s 

where one-electron spin-orbital functions with well-defined orbital angular momentum I (as 
well as with projections mi and rus of the orbital angular momentum and the spin onto the 
quantization axis) are used in the set-up of the corresponding CSF basis. Instead of a defined 
total angular momentum j of each single electron, however, here the orbital angular momenta 
and the spins of all electrons are first coupled together separately in order to obtain a total 
orbital angular momentum, L = li + I2 + ..., and a total spin, S = si + S2 + .... The 
symmetry J of the CSF with respect to a rotation of the coordinates are obtained then from 
the coupling J = L + S of the total orbital and spin angular momenta. In LS"— coupling, 
therefore, L and S are (approximately) also good quantum numbers which can be used for 
the classification of the atomic states. And, in fact, such an LS coupling of the atomic levels 
appears to be well realized in Nature, in particular for light atoms and ions with, say, Z ^ 20 
as well as for the low-lying levels of most neutral atoms. 

As pointed out, the two representations (^ and ((2j of the atomic states are almost equally 
suitable to determine the energies and properties of an atom. To derive the LSJ designation 
of an atomic state, for which a representation {a^ } in jj— coupling is given, it is sufficient 
to expand each CSF \^rJP) into a I/S— coupled basis and to sum up the contributions of the 
different CSF, using the proper weight a^ (t). For two orthonormalized sets of CSF, such a 
transformation 

\-frLSJP) ^^ |7,JP) (3) 

is represented by an orthogonal matrix. The knowledge of this matrix is therefore all that is 
needed in practice in order to carry out the transformation explicitly and to find the symmetry 
of the leading term in the L5— coupled basis. The Lsj program enables the user to evaluate and 
to apply this matrix for any (given) configuration basis in jj— coupling which is used for the 
description of open-shell atoms and ions. There are, however, a few restrictions on the set-up 
of these transformation matrices which are reasonable to maintain for practical computations. 
The next subsection explains these restrictions and summarizes the steps which are needed to 
construct the transformation matrix. 

2.2 Transformation of configuration symmetry functions 

For identical fermions, the symmetry-adapted CSF must represent many-electron functions 
which are also antisymmetric under the interchange of any pair of electron coordinates. In the 



'^Since a jj— coupled basis allows, due to distinctions of equivalent electrons with j — I i: 1/2 (/ > 1), a 
somewhat greater flexibility in the representation of the atomic states, a full equivalence of the representations 
in jj— and LS— coupling is obtained only in the limit of a complete, infinitely large basis. 



atomic (and nuclear) shell model, therefore, it is common practice to construct the CSF from 
antisymmetrized (sub-)shell states which describe N equivalent electrons in a given shell. In 
L5— coupling, for example, such subshell states for the shell (n/) are written as |15j 

\nl^ auLS) , (4) 

where n is the principal quantum number and a represents all those quantum numbers which, 
in addition to the total orbital angular momentum L, total spin S and seniority number i/, 
are needed for the unique classification of these states. In practice, an additional number a is 
required only for shells with orbital angular momenta / > 3, i.e. for electrons from the /-, g~, 
. . . shells. A more detailed discussion of all the possible subshell states (for shells with / < 3) 
has been given recently in Refs. |13| I16j . The CSF with u shells in LS"— coupling takes the 
form 

|7 LS JP) 

= riil^^ n2/2^^..n„/^" aiViLiSi-.-auVuLuSu Li2Si2...Lu~iuSu^iu LS J) . (5) 



The additional intermediate angular momenta Lij, Sij will arise in ^ if the subshell states 
of two or more open shells are coupled to each other as needed for CSF with more than a 
single open shell. 

In jj— coupling, similarly, the subshell states of A^ equivalent electrons of a subshell {uk) are 
represented by ^^ 

= \nK auJ^ (6) 

where k is the relativistic (angular momentum) quantum number 

K = ±{j + 1/2) for I = j± 1/2 (7) 

and where two further quantum numbers v and J are found sufficient in order to classify all 
subshell states with j = 1/2, 3/2, 5/2, 7/2, and 9/2 unambiguously. An additional number 
a is required only for subshell states with angular momenta j > 9/2. Again, Refs. |131 I16j 
provide a more thorough discussion and a complete list of all allowed states for (relativistic) 
subshells with j < 9/2. Of course, the transformation of the subshell states lies at the heart 
of any transformation which is to be carried out for the symmetry-adapted CSF and which 
need to be solved first. The jj <-> LS transformation of the subshell states (in both directions) 
was therefore the topic of a recent extension of the Racah package |141 I17j which was utilized 
also for the implementation of the present program. In going from a jj— to a L5— coupled 
basis, however, there are two difficulties in the transformation procedures which one must 
keep in mind. These difficulties are related (i) to the group-theoretical differences of the 
subshell states in jj— and L5— coupling, respectively, as well as (ii) to the re-coupling of the 
antisymmetrized subshell states which is generally required in order to bring these states into 
a pre-defined coupling order on both sides of the transformation matrix elements. We will 
briefly discuss, in turn, these two issues to make the reader aware also of the limitations of the 
present program. 

Although, in LS— and jj— coupling, the basic notations Q) and © of the subshell states 
appear to be rather similar, these states generally belong to different irreducible representations 
of the 5*03 rotation group because, in jj— coupling, each (nonrelativistic) n/-shell can be 'split' 
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into two (relativistic) subshells with total angular momenta j± = / it 1/2. In the expansion 
of a single shell state \nl^aLS) into a jj— coupled basis, therefore, both subshell states 



with j_ and j-(_, i.e. 




and 




will occur. Here, again, we make use of the 



relativistic quantum number k to simplify the notation below. Obviously, N = N + N and 
K= —{k +1) > must apply where the notation k > and k < has been chosen to be 
consistent with the corresponding total j± = / it 1/2 and subshell angular momenta. 

Making use of this notation, the transformation between the jj— and LS— coupled subshell 
states takes the form 



+ 



(k uJ, k vJ) J 



E 



z(^+^)qLSJ 



- + N 

l(N+N) 



+ 



aLS J I (k vJ, k vJ) J 



(8) 



including a summation over a as well as the total subshell orbital angular momentum L and 



the subshell spin S. As is clear from our notation in Eq. 



(k vJ, k vJ) J 



represents a 



coupled state with well-defined total angular momentum J which is built from the correspond- 
ing jj— coupled subshell states with j± = I ^ 1/2 and the total subshell angular momenta J 

+ 
and J, respectively. This expansion of the jj— coupled subshell states therefore shows that, for 

any number of open (sub-)shells in the construction of the CSF, the transformation of these 
symmetry functions can always be traced back to the transformation coefficients 



_N _- 



+ 



1^^+^^ aLS J\ {k i^J, k I J) J 



(9) 



of the individual subshell. And, in fact, these (Fourier) coefficients are the 'building blocks' 
for any transformation of the atomic states into a L^*- coupled basis. They are implemented 
explicitly into the Lsj component (by means of a recent revision of the Racah package J14j ') 
since their computation requires a recursive — and rather expansive — procedure which traces 
back the subshell states on both sides of expression Q to their parents states with A^ — 1 
electrons, and so on. The transformation coefficients © are implemented for all (sub-) shells 
with / < 3 and occupation numbers A^ = 1, 2, ..., 21 + 1; for all higher (allowed) occupation 
numbers A^ = 2/-|-2, ..., 41+2, these coefficients are simply obtained by using the electron-hole 
symmetry and the tabulated coefficients for A^' = 4/ + 2 — A^ [Tl] . 

A second difficulty in the jj — LS transformation of atomic states arises from the (re-) coupling 
of the antisymmetrized subshell states. Since, in general, CSF with several open shells occur 
very frequently, a pre-defined coupling order of the individual subshell states nili ^ , 712^2 ^ , ... 
need first to be specified, such as 



,N, 



,N2 



N3, 



{■■■{{{li ^aiLiSi, I2 '^a2L2S2)Li2Si2, 1^ ^a3L3S'3)Li235'i23)---) J 



(10) 



for an LS"- coupled basis. Of course, a similar coupling order of the subshell states could 



be written also for a jj— coupled basis 



.A^i 



.N2 



1^2 J 2 



where any (predefined) 



"1 ^1 Ji 

sequence would give rise to a valid many-particle basis. But the predefined sequence of 
jj— coupled subshell states also decides how many (re-) coupling steps are needed to bring 
both sides in the overall transformation matrix (7^ LS JP \ ^rJP) into the same coupling 



order. For this reason, we request the jj— coupled CSF in the Lsj program to be given in a 
(so-called) standard order 

- + ~ + - + 

_JVl_ - +^1+ + _JV2_ - , +^2+ + _JV3_ - , +^3+ + , 

(...(((((ki i/lJl, Ki UiJi)Ji, K.2 I^2J2)Ji2, 1^2 V2J2)Jl2K,z '^3J3)Jl23, 1^3 l^3J3)Jl23---)J ) (H) 



in which these functions have to fulfill two additional conditions: 

(i) If both subshells with common li, i.e. kj and Kj, appear in the expansion, these two subshells 

+ 

always occur successively in the sequence (kj z^jJj, Kj i^iJi) Ji- Formally, we can use this 
sequence even for subshell states with zero occupation if we interpret \k^v = J = 0\ = 1; 
in this case, the full Clebsch-Gordan expansion for the symbolic expression ()lip remains valid 
due to the orthonormality properties of the Clebsch-Gordan coefficients. 

(ii) For the jj — LS transformation of configuration states 
/ (...((Gf^aiLiSi, l^'a2L2S2)Li2Si2. l^'a^L^S^)Li2^Si2^)-) J\ {...{{{{C^T^iJi, ^^'iiJi)Ji, 

_JV2_ - ^ +JV2_,_ + _N3_ - ^ ^JV3+ + \ 

K2 V2J2)J'i2j 1^2 V2J2)Jl2 1^3 VzJ'i)Jl21,i I^Z V3J3)Jl23---) J ) (12) 

we further assume in standard order that li =ki= —{ki +1), I2 =K2= —{k2 +1), •.., i.e. 
that the sequence of (sub-)shell states is the same on both sides of the transformation matrix. 

For such a standard order, explicit expression for the jj — LS transformation matrix of the 
CSF with two open shells in LS— or up to four open shells in jj— coupling are displayed in 
Ref. ^1|. With the present Lsj program, we now extend these expressions to allow for any 
number of open subshells. 

Although, of course, the standard order sets a restriction onto the set-up of the jj— coupled 
CSF basis, it is the order which is typically applied in all relativistic structure codes such as 
Grasp92 PQ or Ratip 0. But even if a non-standard order has been used, the Lsj component 
can still be utilized with only minor additional effort. Since, in fact, the coupling of the subshell 
states is independent of the radial wave functions, it is usually sufficient to (re-) generate 
the CSF basis in standard order and just to repeat the diagonalization of the Hamiltonian 
matrix. Of course, the representation {or (t)} of the atomic states in (^ must always be 
consistent with the basis itself. Moreover, if in the CSF basis the occupied shells appear in a 
different sequence than the usual (non-relativistic) one Is, 2s, 2p, 3s, ..., the condition of the 
standard order is nevertheless fulfilled and the Lsj program can be applied without further 
modifications. In this case, however, the order of the parent states may become different and 
may finally result in a different LSJ classification, if the parent states are accounted for in the 
overall designation of the atomic levels. 

3 Program organization and distribution 

Although Lsj is provided as an independent program for the transformation of atomic levels 
and wave functions from a jj— coupled into a L5— coupled basis, it is designed also as a new 
program component of the Ratip package j^Ell^ and, hence, follows design principles similar 
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to the other components of this package. During recent years, Ratip has been found to be 
a powerful tool for studying the (relativistic atomic) transition and ionization properties of 
atoms and highly-charged ions within rather arbitrary configurations. Apart from the transi- 
tion probabilities and lifetimes, for incorporating the interaction of atoms with the radiation 
field, Ratip now includes further components which help calculate Auger and photoioniza- 
tion properties as well as the coherence transfer in atomic cascades. As described previously, 
the Ratip package is built on Grasp92 ^ which is used to generate the (relativistic) wave 
functions within the MCDF model. 

However, here we shall not recall much about Ratip itself, neither about its overall structure 
nor about previous applications, since this has been described previously 1 18,. Apart from a 
few further applications to incorporate atomic processes with one electron in the continuum, 
our recent concerns pertained to a long life-cycle of the code, making use of an object-oriented 
design and, in particular, the concepts of Fortran 90/95. In fact, all components of Ratip are 
now well prepared to deal with large wave function expansions (^ as needed for the accurate 
handling of open-shell systems. 

3.1 Definition of derived data types 

The proper choice of derived data types certainly belongs to the key features for a long-term 
maintenance and the development of a code. In Ratip, therefore, we introduced derived types 
at various levels of complexity, including for example a data structure cs_function for the 
representation of a single CSF (in jj— coupling), or the derived type cs_basis in order to 
'store' the information about a full configuration basis. To these data structures, memory 
is allocated only on the basis of the actual requirements concerning the number of occupied 
shells and CSF in the given CSF basis. In Lsj, we now followed very similar lines and define 
the corresponding data structures also in L5— coupling. Two simple examples are 



type, public :: cs_f unction_LS 








integer (kind=i lb) 


: : totalJ 








character (len=l) 


: : parity 








integer (kind=i lb) , 


dimension ( 




pointer : 


: occupation 


integer (kind=i lb) , 


dimensionC 




pointer : 


: seniority 


integer (kind=i lb) , 


dimensionC 




pointer : 


: w 


integer (kind=i lb) , 


dimensionC 




pointer : 


: shellL 


integer (kind=i lb) , 


dimensionC 




pointer : 


: shells 


integer (kind=i lb) , 


dimensionC 




pointer : 


: shellLX 


integer (kind=i lb) , 


dimensionC 




pointer : 


: shellSX 



end type cs_function_LS 

to represent a LS"— coupled CSF in terms of the occupation and angular momentum quantum 
numbers of all individual subshells as well as in terms of the intermediate angular momenta 
-^I2)5'i2, -^123) 'S'i23, ... which are stored in the two arrays shellLX and shellSX. The data 
relating to all the L5— coupled CSF are then kept together in the derived data type 



type, public 
integer 
integer 
integer 
integer 
type(nl) 



: : csf _basis_LS 

nocsf ! Number of CSFs in the basis, 
nwshells ! Number of Cnonrelativistic) shells, 
nwcore ! Number of Cclosed) core shells, 
number _of_electrons 



dimensionC:), pointer 
typeCcs_function_LS) , dimensionC:), pointer 
typeCparent_f rom_j j) , dimensionC:), pointer 
end type csf_basis_LS 



shell 

csf 

parent 
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■| lsj_initializatioii_LS_j j I 
I 1 

I 1 

■| lsj_control_transf ormationO I 



data input 



formation of all CSF in LSJ coupling 



jj-LS transformation and classification 
of all selected levels 



data output 



Figure 1: Schematic flow cliart of tlie Lsj program. 



which, within the pointer array parent, also includes details about the parentage of the 
corresponding configuration states functions in the original, jj— coupled CSF basis. Fur- 
ther derived structures concern the representation of the atomic states (asf _function_LS and 
asf _basis_LS) as well as a list of all LS"— coupled subshell states subshell_term_LS to iden- 
tify the allowed terms for open-shell configurations. Obviously, the efficient use of such data 
types facilitates the access to and the handling of selected data. In the Lsj program, all new 
derived types are defined in the module rabs_lsj (file rabs_lsj .f90) which also contains the 
transformation coefficients Q for the jj <^ LS transformation of the subshell states. These 
matrix elements are kept in the derived structure 



type, private ;: LS_jj_ine 

sequence 

integer (kind=i lb) 

integer (kind=i lb) 

integer*8 
end type LS_jj_nie 



w, Q, L, S, J, Nm, Qm, Jm, Qp, Jp 

factor 

nom , denom 



within a integer representation, from which the value of the matrix element follows by com- 
puting simply: fact or y^^. 

3.2 Interactive control and output of the program 

Like the other components of Ratip, the Lsj program is controlled interactively. Figure 
n displays the basic steps in which the program is executed. Following the initialization of 
several data arrays at the beginning of the execution, the control is then taken by the procedure 
lsj_control_transf ormationO, which reads and interprets all input data, carries out the 
transformation for all selected levels, and finally prints the results to screen. This procedure 
also prints a full expansion of all jj— coupled CSF in the corresponding LS"— basis in the file 
input-csl.LS. 
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LSJ: Transformation of atomic states from a jj -coupled into a LS-coupled CSF basis (Fortran 95 version) 
(C) Copyright by G Gaigalas, T Zalandauskas and S Fritzsche, Kassel and Vilnius (2003). 

Transform one or several ASF from a GRASP92 calculation into a LS-coupled CSF basis. 
The transformation starts from the given .csl and .mix files and is carried out 
for the n leading CSF in the jj-coupled basis; the new representation in the LS basis 
is printed similar as in a standard GRASP92 computation. 

Enter the name of the GRASP92 configuration symmetry list file: 
> 
Loading configuration symmetry list file . . . 
There are ... relativistic subshells; 
there are ... relativistic CSFs; 
. . . load complete. 
Enter the name of corresponding .mix mixing coefficient file: 
> 
Loading mixing coefficients file . . . 

. . . load complete; 
Maximum number of considered ASF is: 

Enter the level numbers of the ASF which are to be transformed, 
e.g. 13 4 7-20 48 69 - 85 : 
> 

Enter the number of leading CSF to be printed: 
> 

Figure 2: Interactive dialog for Lsj at the beginning of the transformation. 

Figure El shows the interactive dialog of Lsj which requests the user to provide the Grasp92 
configuration symmetry list ( . csl) file and the wave function representation (of the considered 
level) in a corresponding mixing coefficient (.mix) file. However, independent of the number 
of ASF which are represented in this . mix file, the user then needs to specify the level numbers 
for which the jj — > LS transformation is to be carried out, as well as the number of leading 
L5— coupled CSF which are included in the printout. As discussed before, only the leading 
CSF is used for the LSJ spectroscopic designation of the level but the 'weights' of the next 
most important configuration states may still provide insight into the configuration mixing for 
the given atom or ion. 

As the transformation of atomic levels appears mainly as a geometrical task, the definition of 
the configuration basis along with the corresponding mixing coefficients provide together all 
necessary information for a LSJ classification of the atomic levels. For the transformation 
of the corresponding wave functions, in contrast, some care must be taken also for the radial 
part of the wave function representation, for instance, by taking the average for each pair of 
radial orbitals with angular momentum quantum numbers n > and — k — 1, respectively, 
for all occupied shells (n/). Using the Lsj component, the transformation of the atomic 
levels is carried out promptly for most open-shell structures; the results are therefore printed 
immediately to screen. Moreover, to provide a test of having obtained a correct and complete 
transformation, the sum of all 'weights' in the LS*— coupled basis, i.e. ^^ |Cj (t)|^ (which, 
of course, should be 1), is also printed. 

3.3 Distribution and installation of the program 

Even though the main effort in developing the Lsj component certainly concerned the two 
(new) modules rabs_lsj and rabs_lsj_data, we re-distribute the Ratip package as a whole, 
owing to the tight binding of the new components with the previously developed code. Overall, 
therefore, the present version of Ratip now contains 18 modules as well as the source code for 
the executables of the seven components Anco, Cesd99, Lsj, Rcfp, Relci, Reos99, and 
Utilities. As previously, they are all contained in the ratip root directory. This root also 



13 



Table I. Leading CSFs and corresponding weight coefficients (|cr (r)p and 
|ci ('i^)P) of the low— lying levels of the ls^2s2p configurations in beryllium. 



|*.(^^)) 



n 



weights 



0.66673 



0.33327 




|((2.i)i,(2p-i)i)l) 
|((2.i)i,(2pi)|)l) 



Soo) 



EJc,^''^WP = l 



0.66673 
0.33327 



|((2.i)i,(2pi)|)l) 
|((2.i)i,(2p-i)i)l) 




|*r(J^)) 



L5 



CSFs 



|((2si)2^,(2pi)2p) 3p^) 



,{LS) 



Y.^\cr'ir)? = l 




((2si)2^,(2pi)2p) iPi) 



comprises several makefiles for generating the executable as well as a number of test suites for 
a few components. From this root, the executables of the different components are generated 
by running the command make -f make-component, that is make -f make-Is j in the present 
case. 

By utilizing a large set of different makefiles, it has been found convenient to declare several 
global variables for the compilation and linkage of the program. These declarations are now 
comprised within the script file make-environment which avoids the user having to adopt each 
makefile independently to the local architecture and computational environment. In fact, the 
script make-environment only contains a very few lines but must be sourced before, i.e. by 
typing source make-environment, to make this information available to the system. Using 
this set-up of the code, Ratip has been found portable rather easily to different platforms 
such as IBM RS/6000, Sun OS, or to the PC world. Moreover, a Read. me file in the ratip 
root directory contains further details to the installation. 

For most components of the Ratip package, we provide in addition a test suite in a subdirectory 
test-component of the root where component refers to the names above. For example, the 
directory test-lsj comprises all four (input and output) files in order to run and compare 
the two examples below. These tests are typically designed to exhibit the capabilities of the 
program but may serve also as a quick test of the local installation. 

4 Examples 

To explain the use and execution of the LSJ component, we have selected two examples which 
concern (i) the well-known splitting of the beryllium terms into a singlet and triplet system 
and (ii) shows how the program can be used for the classification of the lanthanides with an 
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Enter the name of the GRASP92 configuration symmetry list file: 

>be-csl.inp 
Loading configuration symmetry list file . . . 
There are 4 relativistic subshells; 
there are 4 relativistic CSFs; 
. . . load complete. 
Enter the name of corresponding .mix mixing coefficient file: 
>bc-rclci.rmx 
Loading mixing coefficients file . . . 

. . . load complete; 
Maximum number of considered ASF is: 4 

Enter the level numbers of the ASF which are to be transformed, 
e.g. 13 4 7-20 48 69 - 85 : 
>2 4 

Enter the number of leading CSF to be printed: 
>2 
Weights of major contributors to ASF in jj-coupling: 
Level J Parity CSF contributions 

2 1 - 0.66673 of 2 0.33327 of 3 

Definition of leading CSF: 

2) 2s ( 1) 2p-( 1) 

1/2 1/2 
1 
Total sum of weights is: 1.000000000000000 
The new level is under the investigation. 
Weights of major contributors to ASF in jj-coupling: 
Level J Parity CSF contributions 

4 1 - 0.66673 of 3 0.33327 of 2 

Definition of leading CSF: 

3) 2s (1) 2p ( 1) 

1/2 3/2 
1 
Total sum of weights is: 1.000000000000000 

Weights of major contributors to ASF in LS-coupling: 
Level J Parity CSF contributions 
2 1 - 1.00000 of 3 

Total sum of weights is: 0.9999999999999996 
4 1 - 1.00000 of 2 

Total sum of weights is: 0.9999999999999991 

Definition of leading CSF: 

2) 2s( 1) 2p( 1) 

2S1 2P1 IP 1 

3) 2s( 1) 2p( 1) 

2S1 2P1 3P 1 
LSJ complete ... . 

Figure 3: Dialog and output for determining the LSJ notation of fow-lying levels of the ls^2s2p 
configurations in beryllium. 
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open /—shell. In both cases, the configuration symmetry list and the corresponding mixing 
coefficients have been generated by means of Grasp92 pQ and the Relci component :4] of 
Ratip, using simply a configuration interaction approach. The input and (expected) output 
of these two examples are provided also in the subdirectory test-Is j. In both cases, the 
transformation is performed very rapidly within less than a second of CPU time on a Pentium 
166 MHz processor. 

For beryllium, the low-lying excited spectrum of the ls'^2s2p configuration consists of just 
four levels which belong to the two ^P and ^P terms and have total angular momenta J = 1 
and J = 0,1,2, respectively. Since the Hamiltonian matrix is diagonal in the total angular 
momentum J, only the two levels with J = 1 will exhibit a configuration mixing and, in fact, 
this interaction is necessary to represent the two levels in terms of the ls'^2s^2p-^ and ls'^2s^2p^ 
relativistic configurations. For the two levels with J = and J = 2, the assignment of the LSJ 
notation is trivial because there is each just one CSF with a weight coefficients equal to 1 in the 
expansion. But to specify also the LSJ symmetry of the remaining two J = 1 levels, we can 
run the Lsj program as shown in Figure Ol where use is made of the configuration symmetry 
list be-csl.inp and the corresponding mixing coefficient file be-relci .mix. As seen from 
the Figure 01 and Table I where the results of the transformation are summarized, the lower 
one of the J = 1 levels (No. 2) has a ^P symmetry, and the ^P term is assigned to level 
No. 4. Moreover, while a strong mixing of the configuration states obviously appears in the 
jj— coupled basis, a (nearly) pure LSJ symmetry is found for all four levels in a L5— coupled 
configuration basis. 

For simple shell structures, such as the low-lying excitation spectrum of beryllium and for many 
other light elements, the LSJ assignment can often be made just on the basis of the excitation 
energies (with respect to the ground state) and total angular momenta of the individual levels, 
if these data are compared with experiment. However, such a procedure is not feasible if two 
or more open shells occur in some electron configuration, or if open d— and /—shell elements 
are considered. In such cases, several tens or even hundreds of levels may belong to a single 
(nonrelativistic) configuration, making a LSJ symmetry assignment on the grounds of the 
calculated energies very unreliable. As a typically example for such a complex configuration, we 
consider the ground state of gadolinium (Z = 64), which has a [Xe] 4/'^5d6s^ configuration with 
3106 possible (jj— coupled) CSF. To determine the LS symmetries of the computed ground 
state with its known total angular momentum 2~, we again run the Lsj component using the 
two files gd-csl.inp and gd-relci.mix, respectively. Figure E] displays the corresponding 
dialog and demonstrates that the ground state of gadolinium has a ^D2 symmetry. The 
leading LS- and jj- coupled CSFs together with their weights (squares of mixing coefficients) 
are presented in Table II. The rather small mixing of other LS terms shows, moreover, that 
the low-lying levels of most (neutral) medium and heavy elements are still well classified in 
terms of their LSJ spectroscopy notation. Despite this classification, however, much larger 
(or even huge) wave function expansions are typically required in order to obtain accurate 
predictions of the level structure and properties of these systems. 

5 Conclusions and outlook 

The new component Lsj of the Ratip package is presented for deriving the LSJ designation 
of atomic levels and wave functions from relativistic calculations. Such spectroscopic notation 
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Enter the name of the GRASP92 configuration symmetry list file: 

>gd-csl.inp 
Loading configuration symmetry list file . . . 
There are 22 relativistic subshells; 
there are 3106 relativistic CSFs; 
>gd-relci.mix 

. . . load complete. 
Enter the name of corresponding .mix mixing coefficient file: 

Loading mixing coefficients file . . . 

. . . load complete; 
Maximum number of considered ASF is: 11 

Enter the level numbers of the ASF which are to be transformed, 
e.g. 13 4 7-20 48 69 - 85 : 
>1 

Enter the number of leading CSF to be printed: 
>30 

Weights of major contributors to ASF in jj -coupling: 
Level J Parity CSF contributions 

1 2 - 0.23193 of 1403 0.21537 of 457 0.12535 of 2443 



0.09234 of 57 0.05331 of 987 



Definition of leading CSF: 

1403) 4f-( 3) 4f (4) 5s ( 2) 5p-( 2) 5p ( 4) 5d-( 1) 
9/2 8 3/2 

7/2 7/2 7/2 7/2 
Total sum of weights is: 0.9999999993812160 

Weights of major contributors to ASF in LS-coupling: 



3s ( 2) 



Level 


J Parity 


CSF contribut 


ions 










1 


2 


0.97616 of 


2 


0.01538 


of 


4 


0.00483 of 


6 






0.00297 of 


1 


0.00034 


of 


8 


0.00005 of 


100 






0.00005 of 


60 


0.00004 


of 


98 


0.00004 of 


58 






0.00003 of 


3 


0.00003 


of 


5 


0.00002 of 


102 






0.00002 of 


62 


0.00001 


of 


15 


0.00001 of 


28 



Total sum over weight is: 0.9999999993812226 



Definition of leading CSF: 

1) 4f( 7) 5s( 2) 





8S0 


ISO 


2) 


4f( 7) 


5s( 2) 




8S0 


ISO 


3) 


4f( 7) 


5s( 2) 




6P0 


ISO 


4) 


4f( 7) 


5s( 2) 




6P0 


ISO 


5) 


4f( 7) 


5s( 2) 




6P0 


ISO 


6) 


4f( 7) 


5s( 2) 




6P0 


ISO 


8) 


4f( 7) 


5s( 2) 




6P0 


ISO 


15) 


4f( 7) 


5s( 2) 




6D0 


ISO 


28) 


4f( 7) 


5s( 2) 




6G0 


ISO 


58) 


4f( 7) 


5s( 2) 




4D1 


ISO 


60) 


4f( 7) 


5s( 2) 




4D1 


ISO 


62) 


4f( 7) 


5s( 2) 




4D1 


ISO 


98) 


4f( 7) 


5s( 2) 




4D6 


ISO 


100) 


4f( 7) 


5s( 2) 




4D6 


ISO 


102) 


4f( 7) 


5s( 2) 




4D6 


ISO 


LSJ complete 







5p( 6) 5d( 1) 

ISO 2D1 

5p( 6) 5d( 1) 

ISO 2D1 

5p( 6) 5d( 1) 

ISO 2D1 

5p( 6) 5d( 1) 

ISO 2D1 

5p( 6) 5d( 1) 

ISO 2D1 

5p( 6) 5d( 1) 

ISO 2D1 

5p( 6) 5d( 1) 

ISO 2D1 

5p( 6) 5d( 1) 

ISO 2D1 

5p( 6) 5d( 1) 

ISO 2D1 

5p( 6) 5d( 1) 

ISO 2D1 

5p( 6) 5d( 1) 

ISO 2D1 

5p( 6) 5d( 1) 

ISO 2D1 

5p( 6) 5d( 1) 

ISO 2D1 
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ISO 2D1 
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6s ( 2) 
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Figure 4: Dialog and output for determining the LSJ notation for the ground state of gadohnium. 
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Table II. Leading CSFs and corresponding weight coefficients (|cr (r)p 
and \cs (7")P) of the ground state of gadohnium. 



^riJ^),, 


'^riJn),s 


weights 


CSFs 


weights 


CSFs 


\c^'^\r)\^ 




Icf'^WP 
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(((4/-^ 1,4/4 8)i,5d-i 1)2) 
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(4/7 8^,5^1 ^D) 9Z?2) 
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{Af ^P,bS ^D) ^Pa) 
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is needed not only to analyze a great deal of spectroscopic data but also for the interpretation 
of many inner-shell phenomena such as the angular distributions of emitted Auger and photo 
electrons. For 2p — ^ 4s resonantly excited argon atoms, for example, it was shown recently 
J19j that the electron emission can be understood properly only if the decay in the (two-step) 
cascade is treated coherently. However, while the atomic levels in the cascade can be specified 
experimentally rather easily owing to the existence of observed line spectra, program tools 
such as Lsj are needed to understand the theoretical data and to select those levels, which 
have to be involved in the computation. 

The design of Lsj within the framework Ratip enables the user to apply a jj — > LS transfor- 
mation for all atoms and shell structures which can be calculated by means of Grasp92 or any 
other Ratip component. This includes, in particular, open d— and /—shell elements which 
have attracted a lot of recent interest. For atomic fermium with Z = 100, for example, first 
measurements on the low-lying excitation spectrum have recently been carried out in Mainz 
and could be supported by means of extended theoretical computations |2D]. Therefore, we 
hope that the Lsj component might be of use also in improving the available date base on the 
lanthanide {Z = 57, ..., 71) and actinide elements {Z = 89, ..., 103) in the future. 

In the design and set-up of Ratip, our intention from the very beginning was to provide the 
user with an open environment for his or her studies of atomic properties. In the last few 
years, therefore, various requests arose for the further development of the code concerning, for 
example, the radiative and dielectronic recombination of highly-charged ions, electron-impact 
processes or even the study of multi-photon excitation and ionization processes. Very likely, 
however, only a few of these requests can be fulfilled by us in the forthcoming years where our 
attention will first be directed to the radiative capture of electrons, i.e. the reversed process to 
the photoionization. Therefore, in order to extend Ratip along additional lines, any support 
of coworkers or even the independent implementation of new components will be appreciated. 
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